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ABSTRACT 

Qin et al. (2004) have derived a formula of count rates based on a model 
of highly symmetric expanding fireballs, where the Doppler effect is the key 
factor to be concerned. In this paper, we employ the formula to both the Qin 
model and the uniform jet model to study how the rising timescale, Ar^r, 
and the decay timescale, Arg j( j, of a local pulse affect the light curve of GRBs. 
Our analysis shows that they do make contributions to both rising and decay 
portions of the light curve of GRBs. Associated with a local pulse with both 
rising and decay portions, the light curve of GRBs in the rising portion is 
expected to undergo a concave phase and then a convex one, which we call a 
"concave-to-convex" character, whereas that in the decay portion is expected 
to evolve an opposite process which is called a "convex-to-concave" character, 
regardless of being in Qin model or in the uniform jet model. These charac- 
teristics are independent of local pulse shapes and their rest frame radiation 
forms. The diversity of light curves are caused by different forms of local 
pulses and different ratios of Are, r to At^- We study a sample of 86 GRB 
pulses detected by the BATSE instrument on board the Compton Gamma 
Ray Observatory and find that the "concave-to-convex" (expected in the ris- 
ing portion of the light curve) and "convex-to-concave" (expected in the decay 
portion) characters, which result from the curvature effect of spherical surface 
(or the Doppler effect), do exist in the light curve of some GRBs. 
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1 INTRODUCTION 



Light curves of gamma-ray bursts (GRBs) vary drastically in morphology, only less than 10 
percent of the light curve profiles detected by the BATSE instrument on board the Compton 
Gamma Ray Observatory (CGRO) can be categorized as being similar in the overall shape. A 
pulse in the gamma-ray light curve is generally considered to result from an individual shock 
episode. Superposition of many such pulses is believed to create the observed diversity and 
complexity of light curves, which could not be describable by a simple formula. Therefore, 
the temporal characteristics of these pulses hold the key to the understanding of the prompt 
radiation of GRBs. 

Some well-separated pulses are seen to comprise a fast rise and an exponential decay 
(FRED) (Fishman et al. 1994). Many attempts have been made to decompose the com- 
plex light curves into pulses, which can be described by several flexible functions based on 
empirical or semi-empirical relations, and quantify their characteristics (e.g., Norris et al. 
1996; Femimore et al. 1996; Lee et al. 2000a, 2000b; Ryde & Svensson 2000, 2002; Ryde & 
Petrosian 2002; Borgonovo & Ryde 2001; Kocevski et al. 2003, hereafter Paper I). E.g., as 
derived in detail in Paper I, a FRED pulse could be well described by equation (22) or (28), 
and using this model, they revealed some features of the profiles of individual GRB pulses. 

However what physics causes such a observed profiles? According to Ryde & Petrosian 
(2002), the simplest scenario accounting for the observed GRB pulses is to assume an im- 
pulsive heating of the leptons and a subsequent cooling and emission. In this scenario, the 
rising phase of the pulse, which is referred to as the dynamic time, arises from the energizing 
of the shell, while the decay phase reflects the cooling and its timescale. However, in general, 
the cooling time for the relevant parameters is too short to explain the pulse durations and 
the resulting cooling spectra are not consistent with observation (Ghisellini et al. 2000). As 
shown by Ryde & Petrosian (2002), this problem could be solved when the curvature effect 
of the expanding fireball surface is taken into account. 

However, there is no consensus on what effects lie behind the observed pulse shapes. We 
suspect that the pulse shapes might have some relations to their counterparts . Motivated 
by this, we explore in this paper possible explanations for the pulse shapes observed based 
on the curvature effect (or the Doppler effect of fireballs). 

Due to the observed output rate of radiation, the observed gamma-ray pulses are believed 
to be produced in a highly relatively outflow, (e.g., an expanding and collimated fireball), 



based on the large energies and the short timescales involved. In the case of a spherical 
fireball, the Doppler effect over the whole fireball must be at work, and the observed FRED 
structure probably provide the intrinsical information. Qin (2002) has derived the flux func- 
tion based on the model of highly symmetric expanding fireball, where the Doppler effect of 
the expanding fireball surface is the key factor to be concerned, and then with this formula, 
Qin (2003) studied how emission and absorbtion lines are affect by the effect. Recently, Qin 
et al. (2004) presented the formula in terms of count rates. Based on their model, many 
characteristics of profiles of observed gamma-ray burst pulses could be explained. They con- 
vince us that the exponential decay phase of the FRED pulse arises from the time delay of 
different parts of the surface and its rising portion is produced by the width of a local pulse. 
However, some interesting questions, such as what a role the rising and decay portions of 
local pulses play in the observed light curve of gamma-ray bursts, and what are responsible 
for the diversity of the FRED pulses, urge us to investigate in a much detail. 

Jets in GRBs were first suggested by Waxman et al. (1998), and then they were widely 
evoked to explain its spectacular energy release (e.g.,Fruchter et al. 1999 ). Subsequent mul- 
tiwavelength observations of GRBs have been interpreted as evidence for explosions with 
jet-like geometry (Stanek et al. 1999; Harrison et al. 1999). Recently, two phenomenological 
models of GRB jets received widespread attention: a uniform jet model in which the emissiv- 
ity is a constant independent of the angle relative to the jet axis (Frail et al. 2001;Panaitescu 
& Kumar 2001; Rossi et al. 2002; Zhang et al. 2002; Bloom et al. 2003; Zhang et al. 2004a, 
2004b), and a structure jet model in which the emissivity is a function of the angle relative 
to the jet axis (Frail et al. 2001; Rmairez-Ruiz et al. 2002; Lamb et al. 2005). In the case of 
a uniform jet, the formula presented in Qin et al. (2004) could be applicable as long as the 
open angle of the jet is taken into account. 

According to the relativistic fireball model, the emission from a spherically expanding 
shell and a jet would be rather similar to each other as long as we are along the jet's axis and 
the Lorentz factor V of the relativistic motion satisfies < 9j, because the matter does 
not have enough time (in its own rest frame) to expand sideways at such situation (Piran 
1995). 

The detection of very energetic photons (in the GeV range) in several gamma-ray bursts 
(GRBs) and the very short variability timescale (Fishman et al. 1995), sometimes less than 
1 ms, exhibited by the 100 keV emission of many bursts have led to the conclusion that they 
arise from sources that are moving at extremely relativistic speeds, with very large Lorentz 



factors T that could exceed 100 (Fenimore, Epstein, &Ho 1993; Piran 1999). At the same 
time, From the derived parameter, 9j, of the sample of 41 source in Friedman et al. (2004), 
we find that the range of 9j is {0.03107, 0.57072 } rad, and its average value is 0.15925 rad. 
The data show that < 9j if T > 100. In this paper, we assume that the condition of 
r _1 < 9j is satisfied for all the source concerned. 

This paper is organized as follows. In section 2, we present the characteristics of the 
light curve of GRBs obtained from our theoretical analysis based on a uniform jet and a 
spherical fireball model. In section 3, we examine the characteristics of light curves of a 
group of source detected by the BATSE instrument on board the Compton Gamma Ray 
Observatory. Discussion and conclusions will be presented in the last section. 



C ( T ) = i„r>2T^i ^2/1 I ' W 



2 NUMERICAL ANALYSIS 

As derived in detail in Qin et al. (2004 ), the expected count rate of the fireball within 
frequency interval [u\, v 2 ] can be calculated with equation (21), which is as follows: 

2nR S jre /(re)(1 + _ r + Te)dT& ^ 9 _^ Adv 
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hcD 2 T 3 (l- (3) 2 (l + kry 
where I(r$) represents the development of the intensity magnitude in the observer frame, 
called as a local pulse function, go,v(v Q 9 ) describes the rest frame radiation mechanisms, and 
r is confined by 1 - cos 6> min + (1 - /5 cos 9 miQ )T 0Mll < r < 1 - cos 9 max + (1-/3 cos 6» max )re imax , 
while f eMn and ?e iinax are determined by f 6Mn = max{r 9imin , T^p^i^ ) and W = 
min{re im ax, T^^i^}^ respectively, we will employ equation (1) in the following analysis. 

2.1 In the case of radiation emitted from the whole fireball surface 

If the radiation is emitted from the whole fireball surface, we can take 9 m i n = and 9 max = 
7r/2. Therefore the range of r is (1 - (3)T e>min < r < 1 + T 6jmax , while Tg Mn and f e>max are 
determined by f fl>min = max{T dMn , r - 1} and r ejmax = mm{r ejmax , r/ (1 - /?)}, respectively. 

Formula (1) suggests that, except the state of the fireball ( i.e., T, R c and D), light curves 
of sources depend on I(t$) and go :U (^o,e)- We assume in this paper the common empirical 
radiation form of GRBs as the rest frame radiation form, the so-called Band function (Band 
et al. 1993) that was frequently, and rather successfully employed to fit the spectra of 
sources (see, e.g., Schaefer et al. 1994; Ford et al. 1995; Preece et al. 1998, 2000). We focus 
our attention to the roles that the rising and decay portions of local pulses play in the light 



curve of gamma-ray bursts. Several forms of local pulses would be employed to study this 
issue. 

(a) Local pulse with a power law rise and a power law decay 
A power law form of I(tq) is adopted as follows: 

I(r e )=I { ^ *- W> (2) 

(r,o < r e < r, max ) 

For this local pulse we find that Arg^wHM = (1 — 2~ 1 ^)(T<j )max — re jm i n ) and therefore, 

_ ^Tg,FWHM _|_ 
'8,max (1— 2~ 1 /m) '0,min- 

For a given value of Tg tmax , the rising timescale, Atq^, and the decay timescale Arg^ of a 
local pulse vary with the shift of the position of r^o, where Arg >r and Ar^ are the FWHM 
widthes in the rising and decay portions, respectively. We assign 

re,o = 70,min + e{r^ max - T e , min ) , whereO < e < 1. (3) 
We find that Ar e>r = r gfi - Te, min , Ar M = T e ^ max - Tg >0 . Let 

r rd {e) = -r- 2 - (4) 
Ar^ 

Quantity r r( i(e) denotes the ratio of the rise timescale to the decay timescale of a local pulse. 

To describe different contributions of Ar r and Atq^ of a local pulse to the light curve, 
we define 

2kR* j£— 7(r fl )(l + /3t,) 2 (1 - r + r e )dr, j£ *^di/ 
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where C{t) = C r {r) + Cdij). Quantities C r (r) and Cd(j) denote the contributions of the 
rising timescale and decay timescale of a local pulse to the light curve, respectively. 

Notes that when /i — 1, the local pulse (2) would become a linear rising and a linear 
decay one. 

Let us employ the Band function radiation form with ao = — 1 and (3q = —2.25, to make 
the light curve, within the frequency range of 100< ^/Vo,p<300. In order to obtain a more 
detailed information of a light curve profile, we evaluate the first order derivatives, C' r (r), 



C' d (r), C'(t), and the second order derivatives C"(t), C^'(r), C"(r), of the count rate C r (r), 
C d (r), C(r), respectively. Their profiles are presented in Fig.l. 

As shown in Fig. 1, two phases, A-B and D-E, undergo a concave process as C^ b (t) > 
and C de (t) > 0, whereas other two phases, B-C and C-D, evolve a convex one because of 
C' BC ( T ) <° and c 'cd( t ) <0- We define 



tbc 

to denote the ratio of the two time intervals from concavity to convexity in the rising portion 
of a light curve, where subscripts "c" and "v" represent concavity and convexity respectively. 

We notice that r cv increases linearly with r rd . For example, r cv ~ 0.28 when r rd =0.2, r cv 
~ 0.76 when r rd =0.5, and r cv ~ 1.64 when r rd =l. Similarly, r r increases, from r r < r d to 
r r > r d , with r rd . We notice that there is an interesting point, B, which is an inflexion from 
concavity to convexity in the rising phase of the light curve, corresponding to the peak point 
in the local pulse. Then the rising timescale Ar^ of a local pulse could be obtained by the 
location of point B from the corresponding light curve. 

We find that both the rising timescale and decay portions of the local pulse make a 
contribution to both portions of the light curve. The rising portion of the local pulse produces 
concave phases on both sides of the light curve and creates a spiky profile, whereas the decay 
one yields a convex phase in the rising portion, and a convex-to-concave phase in the decay 
portion of the light curve, and gives birth to a smooth profile. In a word, the combinational 
contributions of both portions of the local pulse lead to a light curve with its rising portion 
evolving from a concave phase to a convex one (the so-called "concave-to-convex" character), 
and with its decay portion undergoing an opposite process (the so-called "convex-to-concave" 
one). 

(b) Local pulse with a Gaussian form 

Here, we replace I(rg) in equation (1) with a Gaussian form, which is as follows: 

I{r e ) = J exp[-(^^) 2 ] (r fl , min < r e ). (10) 

<j 

The profiles of the light curve are presented in Fig. 2. 

We find that a local pulse with a Gaussian form also produces such a profile of light 
curves: a "concave-to-convex" shape in the rising portion and a reverse one in the decay 
portion. The differences are that a Gaussian form makes a smaller curvature than a power 
law one does in the rising portion of light curves either for the concavity or convexity phases, 



and that the value of r cv keeps to be 1.24 because r rd keeps to be 1 for a Gaussian form. 
However the values of r r and r d change with a. 

(c) Local pulse with a 5 function 
Let us consider the case of I{tq) being a local 5 function, which is as follows: 

^( r ^) = !T 5 ( T e ~ r o,o) (re,min < r e < r e>max ) (11) 

tic 

Associated with a local 5 function pulse, the profile of the light curve can be found in 
Fig. 1 of Qin et al. (2004). There is no rising portion in the light curve for neither rising nor 
decay timescale in the S local pulse. But there is a decay portion due to the curvature effect 
of the expanding fireball. 

Besides the three local pulses discussed above, we also study other forms of local pulses, 
such as /i=2, 3 in equation (2), as well as an exponential form of local pulses, and so on. 
The same conclusions hold for all these case (see Fig. 1 and 2). The value of r cv increases 
linearly with r rd at least in the case of these local forms we investigate. The curvature of 
the rising portion also changes with different forms of local pulses and different r rd , which 
are responsible for the diversity of light curves. 



2.2 In the case of radiation emitted from a uniform jet 

If the radiation is emitted from a jet-like geometry with a half-opening angle 9j and a 
spherical surface (a uniform jet), we assume for such a relativistic outflow both the Lorentz 
factor T and the energy per unit solid angle from the jet axis are constant during the GRB 
phase, and assume that the observer direction is the jet axis. At these assumptions, we 
are able to apply equation (1) to the uniform jet taking 9 m i n = and 9 max = 9, where 9 
represents the half-opening angle of the jet. 

With r=10 and 100, 9 min = and 9 max = 1/T, we repeat the same work as does in the § 
2.1, and find the same characteristics of a light curve of GRBs: a "concave-to-convex" phase 
in the rising portion and a reverse phase in the decay portion, except that there exist a 
cutoff tail feature in the decay portion (see Fig.l in Qin et al. 2004). (the cutoff tail feature 
does not appear in a light curve when 9 max is large enough). The characteristics appear to 
be the result of the curvature effect regardless of a spherical fireball or a uniform jet as long 
as its surface, from which the radiation is emitted, is a spherical one. 



3 OBSERVATIONAL DATA 



To test our analysis mentioned above, we examine the profiles of light curves of GRBs de- 
tected by the BATSE instrument on board the CGRO (Compton Gamma Ray Observatory). 
It is already known that pulses of a GRB show a tendency to self-similarity for different en- 
ergy band (see, e.g. Norris et al. 1996). Therefore, we study in this paper only the count 
rate of channel 3. 

For each source, the background is estimated by a fit with a polynomial expression using 
1.024s resolution data that are available from 10 minutes before the trigger to several minutes 
after the burst. The data along with the background fit coefficients can be obtained from the 
CGRO Science Support Center (CGROSSC) at NASA Goddard Space Flight Center through 
its public archives. After subtracting the background counts, we smooth the data with the 
DB3 wavelet with the MATLAB software in the level of the third-class decomposition. We 
select a sample comprising of 86 well-separated pulses from these data, on which almost no 
other small pulses overlap. 

For the sake of comparison, the selected pulses above are normalized and their variables, 
t, are re-scaled so that the peak count rate is located at t=0 and the FWHM point of the 
decay portion is located at t=0.2. To obtain the value of r cv in a pulse, we evaluate the 
first order derivatives, C"(r), of a pulse, and then identify the inflexion from concavity to 
convexity. Applied this method, we receive a list of the value of r cv for our sample pulses, 
which is shown in Table 1. The distribution of r cv is presented in Fig. 2. 

According to Fig. 3, we divide our sample into 3 classes by the range of r cv , which are 
named as class 1 (r cv < 0.5), class 2 (0.5 < r cv < 2) and class 3 (r cv > 2). Their percentages 
to the total number are 5.8, 58.1 and 36.1, respectively. The observed distribution of r CV) 
whose explanation is not clear, may constrain the mechanism of the birth of local pulses in 
both models. 

Light curve and their slopes of three events, (j 07374 (r cv = 0.42), (J 01883 (r cv = 1.69) 
and jj 01733 (r cv = 5.01), are drawn in Fig. 4. One find that the characters of the slopes are 
in agreement with what we obtain in last section. This suggests that, at least for our sample 
sources, the curvature effect (or the Doppler effect of an expanding spherical surface) for a 
uniform jet or a spherical fireball is indeed at work. 



Table 1. A list of the value of r cv for our sample pulses 




Figure 1. Distribution of r cv of all pulses in our sample 
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Figure 2. Plot of the normalized and re-scaled pulses (first row) and their slopes (second row) for several sample source. 



4 DISCUSSION AND CONCLUSIONS 

All analyses in this paper are based on formula (1), which is suitable for describing light 
curves of spherical fireballs or uniform jets. Numerical calculations show that the characters 
of light curves, the "concave-to-convex" one in the rising phase and the "convex-to-concave" 
one in the decay portion, are not affected by the Lorentz factor T and the half-opening angle 

Let us study the impact of rest frame radiation forms on a light curve profile. As a 
general radiation form observed, the Band function, for which, some sets of typical values 
of the indexes would represent certain mechanisms (see Band et al. 1993), will be employed. 
Here we take three sets of values of the indexes: a —-l and /5 =-2.25, a =-l and (3 —-A.5, 
a =-0.5 and /3 =-2.25. The light curves, slopes, r r and r d , are plotted in Fig.5. 

As shown in the figure, different rest frame radiation forms lead to different values of r r 
and ra, but the value of r cv and the profile character mentioned above are not affected. 

Zhang et al. (2002) found that a structured jet has the same temporal evolutions as 
the uniform jet has as long as the relativistic beaming angle 1 /Y is much smaller than the 
viewing angle 9 V , in this situation, the fireball is observed as if it were isotropic. We thus 
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Figure 3. Plot of C(t), C r (r), Cd(r) — T, C"(r) — r, r r , r^-r for the light curve of the power law local pulse determined 
by (2), where a Band function rest frame form is adopted. And we take r=10, fi = 2, Ato^fw h m=2, TQ >m i n =0, Te tmax =2. 
T 0,O= T 0,max*O-5. cto=-i an d /3o=-2.25 ( for first column ), cto=-l and /?o=-4.5 (for second column), e*o=-0.5 and /3o=-2.25 (for 
last column). The symbols are the same as those adopted in Fig.l. 

suspect that the characteristics of light curves in the structure jet may probably be similar 
to those obtained above, which deserves a detailed analysis. 

One may ask whether or not the characteristics of the profile of light curves come from 
the intrinsic behavior. Our analysis shows that it is not the case. Taking /i — 1 in equation 
(2), which become a linear rise and a linear local pulse, the characters are the same as in 
the case of // = 2. The fact shows that the characteristics indeed result from the curvature 
effect, instead of being an intrinsic one. 

We come to the conclusion that, the light curve characteristics, a "concave-to-convex" 
curve in the rising portion and a revers curve in the decay portion, which is independent 
of the shapes of local pulses, exist in the observed light curve of some GRBs (We find 86 
source in our sample). We believe that reverse characters in both portions (Say, a "convex- 
to-concave" curve in the rising portion and a "concave-to-convex" curve in the decay part) 
would never be detected if the radiated surface is a spherical one, regardless of a spherical 
fireball or a uniform jet. The characters could serve as an indicator of the curvature effect 
or the Doppler effect. 

This work was supported by the Special Funds for Major State Basic Research Projects 
("973") and National Natural Science Foundation of China (No. 10273019). 
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